Optimal thickness of rectangular superconducting microtraps for cold atomic gases 
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We study superconducting microtraps with rectangular shapes for cold atomic gases. We present 
a general argument why microtraps open, if brought close to the surface of the superconductor. 
We show that for a given width of the strips there exists an optimal thickness under which the 
closest distance of the microtrap from the superconductor can be achieved. The distance can 
be significantly improved, if the edge enhancement of the supercurrent near edges and corners 
is exploited. We compare numerical calculations with results from conformal mapping and show 
that conformal mapping can often give useful approximate results. 
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Microstructured magnetic traps are useful devices for 
the storage and manipulation of ultracold atoms and 
Bose- Einstein-Condensates [T]. In particular they allow 
studying the interaction of ultracold atoms with solid- 
state surfaces. However, the lifetime of the atomic clouds 
close to the surface is limited by noise radiation from cur- 
rent fluctuations in the conducting surfaces [31 13] . There- 
fore, the use of superconducting materials for such de- 
vices is being investigated [iHTO]. The energy gap of a 
superconductor strongly suppresses noise in the relevant 
frequency range and thus allows significantly longer life- 
times pdUl4] . Also, the use of superconductors promises 
the construction of hybrid quantum systems and coupling 
with superconducting electronic elements like Josephson 
Junctions, SQUIDs, or microwave cavities as a possible 
route for quantum computation [15l [16] . 

One drawback of superconducting microtraps is the 
strong redistribution of current due to the Meissner ef- 
fect. Cano et al have shown that the Meissner effect in 
superconducting microtraps lowers the trap depth and 
decreases the trapping frequencies fT", TT hampering the 
close approach of the atomic cloud to the superconductor 
surface. The closest distance reported so far in a super- 
conducting microtrap has been 20 /.tm |13j . In principle, 
trap depth and trapping frequencies can always be in- 
creased with stronger currents and fields. However, the 
maximum current is limited by the critical current den- 
sity in a superconducting material. As a possible way out 
of this difficulty it has been suggested to utilize the vor- 
tex state of type-II superconductors in which magnetic 
field partially penetrates the superconductor and reduces 
the Meissner eff'ect fTS^, TS] . However, this reintroduces 
noise radiation coming from the normal conducting cores 
of the vortices and mitigates the use of the superconduct- 
ing material |'20]. 

In this work, we study how the Meissner effect lowers 
the magnetic trap depth. The screening current due to 
the Meissner effect becomes particularly strong near cor- 
ners and edges of a superconducting structure [21,. Here, 
we show that this edge enhancement of the screening cur- 
rent in a microtrap formed by a microstrip of rectangular 



cross-section can be utilized to bring a microtrap signifi- 
cantly closer to the surface of the superconductor. This 
allows a stronger atom-surface coupling without the need 
to enter the vortex state. We perform calculations of 
the current and field distributions in the vicinity of the 
microstrip and show that there exists an optimal strip 
thickness at which the closest distance to the surface can 
be achieved. 

In superconducting microtraps cold atoms with mag- 
netic moments are caught in magnetic field minima. A 
magnetic field minimum can typically be created by the 
inhomogeneous field of a conducting wire and a homoge- 
neous external magnetic field perpendicular to the wire 
In such a situation there will always be a mag- 
netic field zero in the plane perpendicular to the wire. 
There is a fundamental reason why magnetic microtraps 
open in the vicinity of the surface of a superconductor, 
which derives from Ampere's law and the Meissner ef- 
fect: the radial trapping frequency is proportional to the 
normal derivative of the tangential magnetic field compo- 
nent near the surface. Due to Ampere's law we have 
1^ = ^^f-- The normal component i3„ is continuous at 
a surface. As S = inside the superconductor, it follows 
= at the surface. Therefore, the radial trapping 
frequency approaches zero when a superconducting sur- 
face is approached and the trap has to open. 

There are two ways out of this difficulty, which will be 
studied in the following. First of all, the above argument 
is true only on lengthscales larger than the London pen- 
etration depth A, because the magnetic field inside the 
superconductor decreases continuously to zero on that 
scale. Secondly, there can be high screening current flow- 
ing in the superconductor, which may create high local 
field gradients. This is particularly true near corners and 
edges. In the following we will study magnetic micro- 
traps generated by a rectangular superconducting (thin 
film) strip and investigate the behavior near the corners. 

The geometry is shown in the inset to Fig. [l] (a) . The 
cross-section of the strip lies in the x-y-plane with a cur- 
rent / flowing through in negative z-direction. A con- 
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stant external magnetic field Boxt is applied within the 
x-y-plane, such that a magnetic field zero is created some- 
where in the x-y-plane outside the strip. In addition, a 
small offset field Bo of 1 Gauss is applied in z-direction 
in order to avoid Majorana spin-fiip losses [1]. The offset 
field turns the magnetic field zero into a magnetic field 
minimum which locally acts as a harmonic oscillator po- 
tential for the atoms. We calculate the current distribu- 
tion in the strip using the numerical method by Brandt 
and Mikitik . In this method the electromagnetic ker- 
nel is discretized on a grid for the strip cross-section and 
inverted in order to obtain the current distribution in- 
side the strip. For the grid we choose an inhomogeneous 
mesh which is more closely spaced near the corners of the 
strip [5T1 [23] . After determination of the current distri- 
bution, the magnetic field distribution is obtained from 
the Biot-Savart law. 

For a given direction of i?oxt we vary the total current 
/ and the magnitude Scxt in order to minimize the dis- 
tance of the magnetic field minimum from the strip under 
the following two constraints: 1) the magnitude of the 
magnetic field at all points on the surface of the super- 
conductor is required to be larger than the trap minimum 
by an amount of lOksT/ iib in order to ensure that the 
atoms cannot escape from the trap. In the calculations 
below we choose T = 200 nK, a typical value that can 
be achieved in cold atom experiments. 2) The current 
density in the strip is required to stay below the criti- 
cal current density jc at all points. The critical current 
density generally depends on the quality of the material. 
Here, we will take jc = 10^ A/cm^ which is a typical 
value for Nb thin films [Mj ■ Using this optimization pro- 
cedure we finally determine the closest distance of the 
trap minimum allowed under the two constraints. 

In order to check our results and study the infiu- 
ence of the penetration depth A we compare our numer- 
ical calculations with the method of conformal mapping 
[25] . Conformal mapping allows calculation of the field 
distribution and the surface current distribution semi- 
analytically for polygonally shaped strips. This method 
is strictly valid only for penetration depth A — > 0. How- 
ever, it gives useful results as long as the trap minimum 
is farther away from the superconductor than the length- 
scale A, as can be seen from the results below. We use 
the conformal map for a rectangular strip that has been 
given by Chen et al [26j . Within the conformal map- 
ping method the current density diverges at the corners, 
which is a consequence of the A — limit. For our cal- 
culations we need to regularize these divergencies. We 
do so by requiring that the surface current density stays 
below Ajc at a distance of order A from the corners as a 
replacement for the second constraint mentioned above. 

In Fig. [T] we show results for a w = 20 fim wide strip. 
The penetration depth was taken A = 100 nm, a value 
appropriate for Nb thin films. Fig.[T](a) shows the closest 
reachable distance d from the strip as a function of film 
thickness h for an external magnetic field applied in x- 
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FIG. 1: (Color online) Closest reachable distance between 
the trap and a strip of width w = 20 fim as a function of film 
thickness h under the constaints discussed in the main text. 
Black solid circles show the results of the numerical calcula- 
tion, red open circles are the results obtained from conformal 
mapping, a) The trap approaches the strip along the j/-axis, 
external magnetic field is applied along the x-axis. Inset: Ge- 
ometry of the strip, b) The trap approaches the strip along 
the X-axis, external magnetic field is applied along the y-axis. 
The inset shows an enlarged region near the minimum. 



direction. In this case the magnetic field minimum always 
lies on the negative y-axis at the position j/Oi i-S- d = 
\yo\ — h/2. The black solid circles show the results of the 
numerical calculation, while the red open circles are the 
results obtained from conformal mapping. Reducing the 
film thickness h the reachable distance initially decreases 
slightly and reaches a minimum at a film thickness of /i = 
1.2 fim with a distance of 7.64 /im from the strip. Further 
reduction of the film thickness leads to a quick increase of 
the closest distance. Both the numerical calculation and 
the conformal mapping method yield similar results. The 
conformal mapping method deviates from the numerical 
calculation on the order of the penetration depth and 
becomes invalid for smaller strip dimensions, as expected. 

Qualitatively, the appearance of an optimum film 
thickness can be understood as follows: when the film 
thickness is large, the current flows mainly at the surface 
of the strip within a layer of thickness A. For a fixed 
current density jc in the strip the total current / then 
varies with the circumference of the strip. Decreasing h 
leads to a slight decrease oi I oc 1 + h/w. If one wants 
to keep the trap minimum at the same position, one has 
to decrease the magnitude of the external magnetic field 
by the same factor 1 + h/w. Then the magnetic field at 



3 




X [/imj y [fim\ 




X [fim] y [MmJ 



FIG. 2: (Color online) Current distribution at optimized 
distance inside a strip of width w = 20 nm and thickness 
h = 2 jj,m. a) The external magnetic field is applied in x- 
direction and the trap minimum lies on the negative y-axis. 
b) The external magnetic field is applied in y-direction and 
the trap minimum lies on the positive x-axis. 



the surface of the strip will increase because the exter- 
nal field has decreased while the field due to the current 
roughly remained constant because of the smaller thick- 
ness. In effect, the trap depth has increased and this 
means that the trap minimum can be brought closer to 
the strip. This explains the cx 1 -I- h/w behavior of the 
minimum trap distance seen at larger values of h. When 
h becomes smaller, the current flow in the strip becomes 
more homogeneous along the y-direction. Then the total 
current decreases more quickly I oc h. In this case the 
trap depth decreases and the trap minimum has to be 
brought farther away from the strip in order to maintain 
the same trap depth. 

The situation changes completely, if the external mag- 
netic field is applied in y-direction and the trap minimum 
approaches the strip from its thin side along the x-axis, 
as shown in Fig. [l] (b). In this situation strong screen- 
ing currents have to flow at the outer edges of the strip 
in order to screen out the external magnetic field from 
the inner part of the strip. Correspondingly, there exist 
much stronger field gradients in the vicinity of the trap 
minimum and the reachable distance initially decreases 
roughly proportional to h. This behavior stops when a 
film thickness of /i ~ A is reached, because then the mag- 
netic field can penetrate the strip and the screening is 
reduced. As Fig.[l](b) shows, the optimum trap distance 
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FIG. 3: Optimum distance as a function of the critical cur- 
rent density jc for a strip with width w = 20 /im and thickness 
h — 0.2 ^m. a) The external magnetic field is applied in x- 
direction and the trap minimum lies on the negative y-axis. 
b) The external magnetic field is applied in y-direction and 
the trap minimum lies on the positive a;-axis. 



is reached at a film thickness oi h = 0.19 fxm with a 
distance of 0.65 ^m from the strip, about a factor of 12 
times smaller than in Fig. [T] (a) . The conformal mapping 
method does not show the increase of the trap distance 
at low values of h. This is due to the fact that in this 
range h ^ X the conformal mapping method is not valid 
anymore. 

In Fig. [2] we show the current distribution inside a strip 
of width w = 20 /im and thickness h = 2 /im under 
optimized conditions. Fig. [2] (a) shows the current dis- 
tribution when the trap minimum approaches the strip 
along the negative y-axis. In this case, the current flows 
mainly along the three edges of the strip away from the 
trap minimum. There is only small current flowing at the 
edge closest to the trap minimum. This is clear, because 
under optimized conditions the magnetic field becomes 
very small close to this edge and only small screening 
current has to flow here. The highest current density 
(jc) is found at the two corners opposite to the trap. 
Fig. [2] (b) shows the current distribution when the trap 
minimum approaches the strip along the a;-axis. In this 
situation, the current flows mainly along the short edge 
of the strip opposite to the trap minimum. The current 
along the long edges of the strip can almost be neglected. 
This is due to the fact that the external magnetic field is 
applied along y-direction, which creates a strong screen- 
ing current at the short outer edges of the strip . The 
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FIG. 4: (Color online) a) Locus of points of optimum dis- 
tance for different directions of the external magnetic field. 
Here, a strip with width w = 20 fim, thickness h = 2 /xm, and 
jc = 3 ■ 10® A/cm^ was taken. Black solid lines (solid circles) 
show the results of the numerical calculation, red dashed lines 
(open circles) are the results obtained from conformal map- 
ping, b) shows an enlarged region near the corner of the strip. 



highest current density is again found at the two cor- 
ners opposite to the trap. The magnetic field gradient in 
the trap center is 4.5 T/m for the case in Fig. [2] (a) and 
9.9 T/m for Fig.[2](b). 

The optimum distance that can be achieved under the 



two constraints above obviously depends on the quality 
of the material via the critical current density jc- In 
Fig. [3] we demonstrate how the optimum distance varies 
with jc for a strip with width w — 20 /im and thickness 
h — 0.2 fiia along the two coordinate axes. For jc — 
3 • 10® A/cm^, for example, one may approach the strip 
4.5 fim from the y-axis and 0.22 fim from the x-axis. 

In practice, an approach from the thin side of the strip 
may not always be feasible. Nevertheless, the edge en- 
hancement of the screening current can still be exploited, 
if the trap is brought close to the corner of the strip. In 
Fig. |4] we show the locus of points of optimum distance 
for different directions of the external magnetic field for 
a strip with width w — 20 /xm and thickness h = 2 /im. 
Both the results from the numerical calculation and con- 
formal mapping are shown. The results differ only by the 
order of the penetration length, as expected. As one can 
see, the optimum distance lies even closer to the strip in 
the vicinity of the corner. This is related to the stronger 
^-1/3 divergence of the surface current in the corners in 
the limit A — > 0, where I is the distance from the corner 



To conclude, we have shown that for a rectangular su- 
perconducting strip of width w there exists an optimal 
film thickness for which a microtrap can be brought clos- 
est to the surface of the superconductor. The closest 
approach can be achieved near the corners or short edges 
of the strip. This is due to the edge enhancement of the 
screening current. Near the corners distances in the sub- 
micrometer range can be realistically achieved. Such dis- 
tances are essential for good control coupling in quantum 
computation ^16) , for coupling of atoms with a supercon- 
ducting qubit via the near-field of a microwave cavity 
mode [TS], or for studying the Casimir-Polder force. 
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Forschungsgemeinschaft via SFB/TR21. 
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